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Energy dependence of 7 *p total cross section is considered. It is known from the HERA data that the cross section grows with 
energy, and the rate of growth is an increasing function of the virtuality Q 2 of the photon 7*. This dependence can be explained in a 
simple and generic way by the Froissart bound. To implement this mechanism quantitatively a theory satisfying s-channel unitarity 
is required. This is achieved by computing the total cross section from the optical theorem, and the forward elastic amplitude from 
the phase shifts in the impact-parameter representation. A recipe to do so in perturbative QCD is discussed, together with an expose 
of the advantage to calculate the elastic amplitude via phase shift, rather than direct sum of all Feynman diagrams. A two-loop 
computation of phase shift and total cross section is presented, and compared with the HERA data with good agreements. 


1 Total Cross Sections 

1.1 Experimental Data 

7 *p total cross section (Jtot{s,Q) measured at HERA0H 
can be fitted by the formula <ro(Q)(s/A 2 ) a ^\ where s is 
the square of the c.m. energy, and Q 2 is the virtuality of 
the photon. The power a(Q) is equal to 0.08 at Q = 0, 
the same as the corresponding power for all hadronic to¬ 
tal cross sections 0 . We shall refer to this experimental 
feature, that the energy dependence of the total cross 
section is independent of the nature of the colliding par¬ 
ticles, as universality. Moreover, a(Q) is an increasing 
function of Q. Asymptotically, total cross sections are 
limited by the Froissart bound not to grow faster than 
In 2 s, so the power-like increase observed must be ap¬ 
proximate, with the effective power of growth decrasing 
to zero at infinite energies. 

1.2 Q-Dependence of a(Q) 

This last remark provides a simple and generic mecha¬ 
nism to explain why a(Q) is an increasing function of Q. 
From universality, we can write cr tot (s, Q) in the form 

vtotfs, Q) = <Jo(Q)f(s'), (1) 

where <tq (Q) depends both on Q and the nature of the 
colliding particles, but /(s'), which governs the energy 
dependence of the total cross sections, is universal and 
independent of the nature of the colliding particles. It 
is a dimensionless function of the dimensionless energy 
s' = s/A 2 (Q). On dimensional grounds, the energy-scale 
parameter A (Q) is expected to increase with Q. 

To satisfy the Froissart bound, the universal function 
/(s') must asymptotically approach In 2 s'. This means 
that at sufficiently high energies, its slope in a In/(s') vs 
Ins' plot must increase with decreasing s'. A larger Q 
leads to a larger A(Q), a smaller s', and hence a larger 
slope a(Q) in the log-log plot. This shows that a(Q) is an 


increasing function of Q, just as observed in the HERA 
data Eli. 

2 QCD Phase Shifts 

2.1 Impact-Parameter Representation 

To be able to implement this mechanism quantitatively, 
we must start from a theory capable of describing rising 
total cross sections which obey the Froissart bound. This 
is so if total cross sections are computed from the optical 
theorem, and the forward elastic scattering ampltudes 
are calculated from phase shifts in the impact-parameter 
representation: 

o’tot(s) = -ImA o (s,0) (2) 

s 

A(s, A) = -is j d 2 be iKX (e 2iS ^ - l) (3) 

The impact parameter b and its conjugate, the momen¬ 
tum transfer A, are two-dimensional vectors in the trans¬ 
verse plane. In QCD, the amplitude A is a colour matrix 
whose matrix elements are the scattering amplitudes of 
partons of various colours. The amplitude Aq appearing 
in (] 2 j) is obtained from A of (||) by retaining only the 
colour-singlet exchange amplitude. 

The effective size R measured in the total cross sec¬ 
tion can be estimated from the condition S(s,R) ~ 1. 
If the interaction has a range /r _1 , so that S(s,b ) ~ 
£(s) exp(— pb) for large b, then the effective radius is given 
by R = /r” 1 ln/(s). This grows with enegy if £(s) does, 
and as long as £(s) does not grow faster than a power, R 
will not rise more rapidly than Ins, and er tot ~ ttR 2 will 
not grow faster than In 2 s, as prescribed by the Froissart 
bound. 

2.2 QCD Phase Shifts 

So far the discussion has been general and qualitative. 
To proceed further we must find a way to calculate the 
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phase shift S(s,b). We shall discuss how to do so within 
the framework of perturbative QCD. 

Scattering amplitudes are given by Feynman dia¬ 
grams. From (|3j), the phase shift is given by the loga¬ 
rithm of the amplitude in the impact-parameter space. 
With the presence of non-commuting colour matrices in 
the amplitude this calculation promises to be compli¬ 
cated. It is not a priori clear that any simple formula 
can be obtained for the phase shift. Fortunately, in the 
leading-log approximation (LLA) where only the high¬ 
est P 0 Wf 1 i'rfi4r'if 1 i s kept, the result turns out to be very 
simple □’BffijoEl We will now state this result for quark- 
quark scattering. 

First, use the SU(N C ) commutation relation 
[A a ,Ab] = ifabc^c for colour matrices A a in any rep¬ 
resentation, together with the normalization condition 
Tr(A a Ab) = 6 a b/ 2 for the colour matrices in the fun¬ 
damental representation, to reduce the colour factor of 
each Feynman diagram into combination of a linearly 
independent set of ‘pseudo-planar’ colour factors, as il¬ 
lustrated in Fig. 1. Choose from this set all the primitive 
colour factors Gi , which are those that remain connected 
when the two quark lines are removed from the diagram. 
Of those shown in Fig. 1, G a and Gd are primitive, 
and Gb,G c ,G e ,Gf are not. Since many Feynman dia¬ 
grams can contribute to the same primitive colour factor 
Gi, the corresponding amplitude di(s,b) in the impact- 
parameter space will generally be given by a sum of many 
Feynman diagrams. 

The phase shift is obtained by summing the contri¬ 
butions from all primitive colour diagrams: 

2 5(s,6) = X)d i (a,6)G < . (4) 

i 

It is clear from |) that a S(s,b) calculate even in 
the lowest order will give rise to an expression for the 
amplitude A(s, A) involving amplitudes of all orders. In 
this way a phase shift calculation of an amplitude is sim¬ 
ilar to a renormalization-group improved calculation of a 
Green’s function or a physical parameter. Moreover, it 
can also be shown B that the g 2 and Ins dependences of 
di{s,b) are of the form g 2m D(g 2 In s), if Gi has m gluon 
lines connecting the two quarks (to = 1,2,3, 2,4,4 re¬ 
spectively for (a), (b), (c), (d), (e) in Fig. 1). Let us 
illustrate the significance of these two remarks. 

Suppose we calculate electron-electron scattering in 
QED, via the exchange of n photons. We will ignore 
all fermion loops in the intermediate states. There are 
n! (crossed and uncrossed) ladder diagrams describing 
this process, whose individual powers of Ins grow with 
n. However, it is known that when the n\ diagrams are 
summed, all powers of In s cancel one other, leaving be¬ 
hind a contribution to the total cross section independent 
of energy. This can be seen from (Q) in the following 


la; (b) (c) (d) (e) (f) 

Figure 1 : Colour diagrams. The colour factors (a) and (d) are 
primitive, the others are not. 

way. In this case, the only primitive ‘colour’ diagram is 
a single-photon exchange diagram, G a , for all the others 
are disconnected once the two electron lines are removed 
(Gd is not present for QED). Among the n! ladder dia¬ 
grams, the only one that would contribute to G a is the 
Born amplitude, with n = 1. Hence the phase shift is 
given by the Born amplitude, which is proportional to g 2 
and contains no Ins. This agrees with the general form 
discussed above with to = 1. Moreover, when we ex¬ 
pand the exponential in (j|) to order g 2n , we see that the 
amplitude to that order must also contain no In s terms, 
agreeing with the result obtained by direct calculations. 

For quark-quark scattering in QCD, gluon self cou¬ 
pling allows an infinite number of primitive colour factors 
to be formed, even when quark loops are omitted. G a is 
colourwise a single gluon exchange, with amplitude of 
the form d a (s,b) ~ g 2 D(g 2 \ns). When all diagrams are 
summed in LLA, this gives rise to the familiar Reggeized 
gluon amplitude with D being an exponential function. 
Otherwise D is a polynomial if the sum is taken only to 
a finite order. When the exponential in ®) is expanded 
with S(s,b) ~ d a G a being the exponent, the constant 
term is cancelled and the linear term has colour-octet 
exchange, so it is starting from the quadratic term that 
it may contribute to the colour-singlet amplitude Aq in 
(j2|). For that reason the contribution from G a to er tot , 
is of order g 4 , times a function of g 2 lns. In the region 
where g 2 is small and g 2 In s is of order unity, the con¬ 
tribution from G a to a to t is then 0(g 4 ). The only other 
terms of the same order came from primitive colour fac¬ 
tors that have two vertical gluon lines between the two 
quarks (to = 2), such as Gd in Fig. 1. Those diagrams are 
also of 0(g 4 ), and since two vertical gluons are involved, 
the primitive colour factors do contribute directly to a 
colour-singlet exchange. The BFKL Pomeron amplitude 
is therefore the 0(g i ) contributions from the m = 1 and 
the to = 2 primitive colour factors. To be sure, other 
terms in the expansions of the exponential give higher 
order contributions, but these subleading-log terms are 
necessary to uphold s-channel unitarity and the Froissart 
bound. 
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3 Two-Loop Phase Shifts 

The quark-quark scattering amplitude in LLA, up to two 
loopa,or 0{g 6 ), can be found in the book of Cheng and 
Wut3. All we need to do is to extract from it the phase 
shift <5(s, b) using eq. (Q), and then the colour-singlet ex¬ 
change amplitude A 0 in (j3). The result_for the universal 
energy function f(s') in (Ty) is given by a 


/(*') =Jd 2 




(26 0 
cos ( 3 


1 —2D ( 

3 V 3 


26 a {s 1 b) = - g 2 K 0 (x ) + (g 4 N c lns'/4ir)K 2 (x) 

-(g 6 N 2 In 2 s'/32n 2 )V(x), 

D(x ) = —i6 d {x) = (g 6 \ns'/4n)[Ko(x) - ^ V(x)\, (5) 

where Kq(x) is the Bessel A'-function, and V(x) = 
-(1/2tt) 4 V2 f d 2 x'K$(\x- x'\)K$(\x'\). 


4 Comparison with the HERA Data 

Equation (^) is computed for quark-quark scattering, but 
if the energy function f(s') is indeed universal, as exper¬ 
imental data indicate, (see the discussion leading to (I])), 
then we can use it to compare with experimental results. 
Even so, the perturbative formula in (|5|) is expected to be 
quantitatively accurate only when both g 2 and g 2 In s' are 
small. The former is so when the virtuality Q 2 is large, 
but the latter is not satisfied for the present data. Higher- 
order calculations are needed and they are underway □, 
but for now, to get a qualitative idea, let us compare 
with data anyhow. Choosing a linear dependence for the 
ennergy scale parameter, A (Q) = Ao + cQ, with Ao = 0.2 
GeV and c = 4, the comparison is shown in Fig. 1. We 
see that inspite of the small Q available in the data, the 
agreement of theory and data seems to be very good. 

This agreement can probably be attributed to two 
things. First, the mechanism explained in section 1.2 is 
so general that any reasonable function f(s') satisfying 
the Froissart bound is expected to yield a reasonable re¬ 
sult, if a suitable energy-scale parameter A (Q) is chosen. 
Moreover, the agreement is reminiscent of the success in 
the early days of QCD, before higher-loop formulas be¬ 
came available. Drell-Yan pair mass spectrum and other 
quantities were correctly predicted, though normaliza¬ 
tions were off by the famous ‘K factor’. When higher- 
order corrections are taken into account, there will be 
a correction to the universal energy function f(s') with 
a corresponding change in the energy-scale parameter 
A (Q) to fit the data. This is not unlike A qcd, which 
takes on different values when deduced from experimen¬ 
tally measured a s (Q ), depending on whether it appears 
in a one-loop or a two-loop formula. 
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Figure 2: 7 *-protop total cross-sections as a function of .s. Data 
are taken from Ref. El where references to the original experiments 
can be found. Q 2 of the photon, in (GeV) 2 , are given on the 
right. Where two numbers are listed, the second one refers to the 
low-energy data of E665. The solid curve is the prediction of the 
two-loop phase shift calculation, eq. (El). The dotted lines are s 0 08 , 
and the dash line is s 0 05 , both there for comparisons. 


5 Conclusion 

Experimentally total cross section grows with energy 
like with an effective power a(Q) increasing with 

the virtuality Q 2 of the photon 7 *. This behaviour can 
be explained by universality and the Froissart bound. To 
compute the universal energy function f(s') governing 
the energy dependence of all total cross sections, we re¬ 
sort to the impact-parameter formalism which preserves 
s-channel unitarity and the Froissart bound. We have 
discussed a recipe in perturbative QCD to calculate the 
phase shift S(s,b) needed in that formalism. We have also 
discussed the advantage of calculating the elastic ampli¬ 
tude using phase shift, rather than doing it directly from 
summing all Feynman diagrams. Formulas are given for 
quark-quark phase shifts to two loops, and for the corre¬ 
sponding energy function f(s'). When they are used to 
compare with the HERA experimental data, good agree¬ 
ments are obtained, the reasons for which are also dis- 
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cussed. 
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